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Abstract 

Instanton analysis is applied to models B — H of critical dynamics. It is shown that 
the static instanton of the massless (p^ model determines the large-order asymptotes 
of the perturbation expansion of these near-equilibrium dynamic models leading to 
factorial growth with the order of perturbation theory. 

Key words: instanton, dynamic models, large orders 
PACS: ll.lO.Jj, 05.70.Jk 



1 Introduction 



Large-order asymptotic behaviour of perturbation series of the paradigmatic 
static (j)^ field-theoretic model has been thoroughly explored with the aid of 
instanton analysis and applied to the resummation of the asymptotic series 
brought about by the standard perturbation expansion [1,2,3]. 

However, large-order asymptotics in dynamic field theories constructed from 
Langevin equations with the use of the Martin-Siggia-Rose (MSR) formal- 
ism [4] is almost completely unexplored. 
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In this paper we use the recently proposed dynamic-instanton method [5] 
to assess large-order asymptotic behaviour of dynamic models from B to H 
in the standard classification [6]. All these models possess Gibbsian static 
limit and thus describe dynamics near thermodynamic equilibrium. In the 
majority of these models mode-coupling effects are present requiring significant 
development of the method used for model A in Ref. [5] 

We construct the stationarity equations for the steepest-descent method and 
give arguments in favor of existence of the dynamic instanton solution. Large- 
order asymptotes of correlation functions, response functions and critical ex- 
ponents are shown to be determined by the instanton solution of the corre- 
sponding static problem, which in all cases considered may be constructed 
with the aid of the known static instanton of the massless model. 

The present article is organized as follows: in Sec. 2 representation of the MSR 
field theory in the functional-integral form for the generic model of critical 
dynamics from the corresponding Langevin equation system is reviewed. Then 
in Sec. 3 we show that dynamic instanton may be constructed and the steepest- 
descent evaluation of the functional integral carried out in the same fashion 
as for model A [5] with the result that the leading contribution is given by the 
corresponding static instanton. Construction of the dynamic instanton with 
the aid of the static instanton of the massless 0^ model for models B — H of 
critical dynamics is analyzed in Sec. 4. Results of this paper are summarized in 
Sec. 5. Calculation of the fiuctuation determinant is reviewed in the Appendix. 



2 Field theory for critical dynamics 

Standard models of critical dynamics may be described by the generic Langevin 
equation 



determining dynamics of fluctuations in a system with the thermodynamic 
equilibrium corresponding to the field theory with the static " action" (effective 
Hamiltonian) S. In (1) the source ^ is a Gaussian random field with zero mean 
and the correlation function 



chosen to maintain the fiuctuation-dissipation theorem. To this end the fol- 
lowing conditions are imposed on the terms of the kinetic coefficient as well: 





(2) 



= 0. 



(3) 
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The symmetric part of the kinetic coefficient a is assumed to be field-independent 
and positive-definite. In equation (f ) the static action is assumed to have three- 
point and four-point interaction terms only and to contain the action of the 
massless model 

S4 = lwf + l^<P\ (4) 

for which the instanton solution is known [1]. Relations (1) - (3) constitute 
a compact form of the description of all standard models of critical dynamics 
[7]. Here and henceforth, all necessary integrals and sums are implied. In the 
forthcoming treatment the full d-dimensional Euclidian space R'^ is assumed 
for spatial variables, but for the time variable it is convenient to carry out most 
of the analysis on a finite time interval [to, T] and pass to the limit to — > — oo 
and T — >• oo at a later stage. 

In the MSR approach, functional integrals for correlation and response func- 
tions are calculated with the "measure" 



Tl[ip, ip'] = VipVip'^Jdet Ml Mo e'^ , (5) 

where the dynamic action for the model (1)- (3) may be written symbolically 

[4] as 



+ [Oiab + Pah, 



(6) 



Ot Oifib 

We stick to the choice of Ref. [5] of the determinant factor in representation 
(5) in the field- independent form with the self-adjoint operator M^Mq, where 
the operator Mq is generated by the free-field part of dynamic action (6) as 

(Mo)„, = + (aac + K^b . (7) 

Here, P^^^ is the field-independent part of the skew-symmetric term of the 
kinetic coefficient, whereas Kab is the kernel of the quadratic form Sq of the 
static action: 

'S'o = 2 Vo.KabVb , (8) 

and is thus, by definition, symmetric and positive. The choice of the determi- 
nant factor in measure (5) implies that the interaction functional Si is taken 
in the normal form, which in the perturbation theory is equivalent to the 
condition that there are no graphs with closed loops of the retarded propa- 
gators [Mq^^^ ^(t,x;t',x') attached to the interaction vertex. This amounts 
to vanishing retarded free-field propagator at coinciding times, and no field- 
dependent factors to the determinant det M^Mo are generated. 

We remind that initially in the ordinary perturbation theory the boundary 
conditions on the time axis fix ip'{T) = 0, but leave ^{T), (p'{to) and (p'{T) 
arbitrary. In the instanton approach the space of integration E{A) in the 
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functional integral is chosen using the properties of the quadratic form of the 
expansion of dynamic action at the saddle point and therefore it is in general 
different from that of the ordinary perturbation theory, the construction being 
similar to that for model A [5] . 



3 Instanton analysis for the generic model of critical dynamics 



We are dealing with a dynamic model which has several coupling constants 
which, however, are proportional to some powers of a single parameter, e.g. e, 
at the fixed point of the renormalization group equation). Therefore, it seems 
most reasonable to base the large-order analysis on the number of loops rather 
than powers of different coupling constants. This may be effected by extracting 
a common factor u from all coupling constants and then counting large orders 
by powers of u. Thus, we cast the interaction part of the static action in the 
form 

Vabc<Pa<Pb^c + gabcd^a^b^c<Pd = UVabc^a^b^c + U^Qabcd^a^b^c^d ■ (9) 

The task then amounts to a steepest-descent calculation of the parametric 
integral 

1 [[v^V^'^,,{t,, xi) . . . Xk) e--^ - ^1"^ 

' , (10) 



e 



'0 



where $ = {99, (/?'}, expressing the A^th order contribution to perturbation 
expansion in u of the /c-point Green function. In representation (10) the nor- 
malizing determinant has been written in the form of a functional integral with 
the free-field part of the dynamic action in order to make explicit cancellation 
of Jacobians related to various changes of variables. As shown in Rcf. [5], an 
instanton solution for expression (10) for model A may be constructed in close 
relation to the instanton solution of the corresponding equilibrium model. In 
this Section we show that a similar procedure allows to reduce the large-order 
analysis of the dynamic functional integral (10) to that of the static problem 
determined by the static action S. 

It is convenient to extract the large factor A^ in the exponential of the numer- 
ator of expression (10) by a suitable change of variables accompanied by the 
same transformation in the denominator to avoid uninteresting Jacobi deter- 
minants: (pa — > ^/N(pa, ^'a V^V'a) ^ ~^ u / ^/N . To simplify the forthcoming 
calculation of the fluctuation integral around the saddle point we also scale 
out the dependence on the variable u hj ipa ^ ipa/{iu), ip'^ — >• ip'J{iu). As a 
result of these transformations we arrive — up to the factor N^^'^^^/'^{iu)~^ 
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at the representation 



27ri / 



where the dynamic action S (as well as the static action S) are independent of 
the loop-counting variable u and these can be obtained from the initial ones 
by change of the coupling constants Qahcd ^ -Qobcd-, ^abc ^ -ivabc- 

The stationarity equations of the method of steepest descent assume the form 

S(fia dt SipcSipa Sipa Slfic 

-— - = -2Q;ab(^ft + — — + (Q!a6 + — = , (13) 

S=4. (14) 

As in static theory, equations (12) - (14) have a nontrivial solution with neg- 
ative u"^ only. The dynamic instanton for the basic field is the nontrivial 
solution of the equation 

^'^^ + {aab - Pab) 1^ = 0. (15) 



Ot d(pi, 

Substitution of this solution to equation (13) leads to the nontrivial instanton 
(fi'j-) for the auxiliary field, whose equation we quote in the two equivalent forms 
used later: 

SS I (d^i SS\ 

With the use of the instanton equation (15) and the antisymmetry property of 
the kinetic coefficient Pat it may be readily seen that the stationarity equation 
(12) is fulfilled on solution (16). 

In addition to equations (12), (13) and (14) variation on the boundaries of the 
time interval produces [assuming fixed initial field </?(io,x)] the condition 

^'(T,x) = 0, (17) 

which - due to the first equation (16) - singles out the solution approaching the 
static instanton at the final time instant. Indeed, according to relations (16) 
and (17) the dynamic instanton is a function conforming to the Cauchy 
condition for the final time instant ipd{T,:k) = ipst{y^), where ipst is the static 
instanton solution for the static action S. It should be noted that the iterative 
solution of the instanton equation (15) also has the property lim (foit, x) = 
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0, because — due to the sign changes in comparison with (13) — the Green's 
function of hnear part of equation (15) vanishes infinitely far in the past. A 
similar equation for the model A was investigated in [5], where it was shown 
that the solution of this equation on the infinite time range for negative 
tends to zero at t — > — oo and to (fst at i — > oo. 

Substitution of the solution (p^ and (p'jj in dynamic action (6) reveals that the 
dynamic action on the dynamic instanton solution gives rise to a time integral 
of a total time derivative (the expression dS/dt = {dip/dt){6S/5ip) was used) 
asymptotically coincides with the static action on the static instanton: 

S i^D, ^'d) = S {iPst) - S {ipo) > S {ipst) , 

to— >— oo 

here (/3o(x) = (poitQ^n). Equation (15) gives rise to iterative solution with 
an advanced Green's function, therefore the field (po vanishes in the limit 
to — > — oo 

Similar the third saddle-point equation (14) on solution (15) and (16) reduces 
to the following equation for the limiting values of the dynamic instanton 

S {(fist) - S {(fio) = -— , 
and thus recovers the stationarity equation for the static instanton 

S{^st) = -Y^ to^-oo. (18) 

Thus, as for model A, we arrive at the conclusion that the exponential factor in 
our steepest descent analysis of Green function (11) as well as pre-exponential 
factor asymptotically are the same as in the corresponding equilibrium static 
theory, provided there is a nontrivial solution of equation (15). Such a solution 
may be constructed in the same fashion as in model A as the usual tree-graph 
solution of the nonlinear equation (15) with the given final Cauchy condition 
ifiD{T,x) = (fist{^) instead of an initial condition. 

In the integration space of the functional integral of the A;-point Green func- 
tions (11) the degeneracies of the instanton due to the usual dilatation in- 

variance and translation invariance in coordinate space may be removed with 
the aid of the same unit decomposition [1,10,11] as in the static instanton 
theory with the corresponding conditions imposed on on the final value of the 
integration field (fi{T, x). 
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Thus, we arrive at the Green function 




1 rdujj I^^W^a, (ti, xi) . . . (4, Xk) I e 



N(^-S/u'^ + \nu) 




(19) 



where I stands for all contributions from the degeneracy-lifting unit decom- 
position. At the leading order in we replace all $(ti,Xi) by {(/?£), in 
the pre-exponential factor. In the stationary limit to ^ —oo the initial value 
of the instanton (pnito) vanishes and in the exponential the static instanton 
action is recovered at leading order in N. 

The most subtle point of the saddle-point approach, i.e. the calculation of the 
fluctuation determinant at the instanton solution, is outlined in the Appendix. 
The main conclusion is ~ as in the case of the model A [5] - that the fluctuation 
determinant coincides with that of the static theory. The analysis of large-order 
contributions to correlation and response functions also proceeds in the same 
way as for model A and we do not dwell on it here. 

Thus, we conclude that in the translation-invariant in time case the asymp- 
totic properties of the dynamic model at leading order in A^ are determined by 
the static instanton solution which leads to factorial growth of the large-order 
contributions as in the static instanton analysis. Therefore, the large-order be- 
haviour of an arbitrary quantity F (correlation or response function or critical 
index) may be expressed as 



where F^^^ is the A^th order contribution to F of the expansion in the param- 
eter e (e is the coupling constant g or the dimensional regularization parame- 
ter e). In (20) the constant om depends on the model only (M = A, B, C, ...), 
whereas the exponent h and the amplitude factor C - either a constant or a 
function of coordinate and time arguments - depend on the quantity F as 
well. See, for instance, the analysis of the factor C for a generic function in 
the static theory presented in Ref. [12]). 

As in the case of the model A [5], when expression (20) is used for contributions 
to correlation or response functions, the $ fields in the pre-exponential factor 

of (19) have to be replaced by the static instanton ipat for the fields ip and 
put equal to zero for the fields ip' at the leading order in 1/A^. The dynamic 
properties of the Green functions are thus determined by the first correction 



F^""^ = C N\a'^N\ 



(20) 



in 1/A^ [5]. 
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4 Static and dynamic instantons for models B — H 



We have shown in Sec. 3 that at the leading order the results of the dynamic- 
instanton calculation are determined by the static instanton corresponding to 
the static action S. In this Section we review the results of the static-instanton 
analysis, which in all cases are based on the long-known static instanton of 
the massless 0^ model. 



4.1 Model B 

In model B only one field is present, the static action is that of the massless 0^ 
model (4) and the symmetric kinetic coefficient is proportional to the Laplace 
operator: 



Since the static action is exactly the same as in model A and the large-order 
behaviour of the perturbation expansion is determined by the static instanton, 
the conclusions about this behaviour are the same as for model A: in model 
B expression (20) is valid with the constant as — ao known from the analysis 
of the static theory [1,10] [for the action (4) = l/(167r^)]. 

Although the explicit form of the dynamic instanton was not used in the anal- 
ysis presented above it would be useful for the investigation of the amplitudes 
C in relation (20) for the dynamic Green functions. Therefore, we remind how 
the dynamic instanton may be constructed. A convergent iterative solution of 
the equation (15) for model B can be written in analogy with model A [5], 
viz. 



can be taken as the zeroth-order approximation. The next orders of the ex- 
pansion are constructed as the usual tree-graph solution of equation (15) with 
the aid of the advanced kernel of the linear operator [dt — AV^), which can be 
represented in the form 



(21) 





(22) 
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4.2 Model C 



In model C there are two components of the basic field: let (pi — (f) and (p2 — m, 
and the static action is of the form: 

sc^l wr + 1*^' + T + ^ ■ ^^^^ 

The kinetic coefficient is still symmetric and diagonal 

(24) 



't ^ 



a 



y -AV2 J 



To keep track of the number of loops we put g — and V2 — xu. Then the 
static instanton is determined by the set of equations 

- - - ixm0 = , (25) 
m-^X(/)2 = 0, (26) 



^ = -1 (27) 



and thus may be constructed with the use of the usual instanton solution ^gt 
of the massless 0^ model: 

-V2</>.t-^0L = O, (28) 



i /cix (V0,,)' 0.* - I /^x = , (29) 



with ^ = - 1. To use large-order asymptotes for investigation of the e 
expansion we write x = x(e) = f 2*/ ( where ^2* and (7* are the coordinates 
of a fixed point of the renormalization-group fiow equations) . 

The stabihty of the fixed point in the model considered depends on the value 

of the critical index of the specific heat a. If the index a < 0, then ^2* = 
at the stable fixed point. In this case g = —\ and the constant a = in 
expression (20) for this model is the same as in the model A and the usual 
static 0^ model [1,2]: ac — Oo- 

In the opposite case, when index a > we have (7* ~ e and ~ e at the 
stable fixed point giving rise to a nontrivial e expansion for the parameter 
X = V2IU [7]: 

x=x(°) + x«e + x(2)e2 + ... (30) 
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Therefore, the constant ac in relation (20) is now determined by the equation 



1 [1 - 3x(e) 



21 



Oo ac 

Taking into account that in the e expansion asymptotic structure (20) appears 
as 

p[N]^N ^ NiCa^N^'e'' (31) 

we see that the only contribution to the constant ac comes from the x*^°^ term 
of expansion (30). In this case, therefore, 

ac = 



1 - 3x(o)2 ■ 



The coefficient x'^^-' from expression (30) determines the correction to the 
amplitude factor C of the form 



exp 



1 - 3x(o)2 



whereas x^^^ and higher-order terms of the e expansion of x(e) are irrelevant 
for the leading order in 1/N. 

Let us note that actually we are considering here a ^/e expansion, since the 
loop-counting parameter u ~ y/e at the leading order, and due to the scaling 
4>,m ^ 4>/{iu),m/ (iu) the u expansion of correlation functions with an odd 
number of m fields contains odd powers of u whereas those with an even 
number of m fields contain even powers of u only. Nevertheless, we may use 
the usual form (31) of the perturbation expansion taking into account simple 
additional normalization in the functional-integral representation (11) for the 
Green functions. 



The fluctuation determinant is also the same as in the massless static model, 
because it is given by the Gaussian integral with the quadratic form 

l^^S^^ 5ct> (- + + ^2m,t) 5ct> + 5ct>V2ct>st5m + ^Sm\ (32) 

Here, the gaussian integral over dm gives the determinant of the unit operator 
and produces the Gaussian integral over 50 with the quadratic form 

Is^S^l 6<p (-V^ + 0,^,) 6<p = \ 50 (- - ^ 0^,) 54> (33) 

which is exactly that of the massless (p^ model. A subsequent integration over 
5u yields the additional factor 1 / Vl — Sx^. 
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Equation (15) for the dynamic instanton here assume the form 

^ = r (- W - - i'^-t') , (34) 

The zeroth-order approximation of the iterative solution can be written as 

(/.g^(i,x) = J dx'G'(T-i,x-x')0.t(x'), 

mg)(t,x) = ^1 rfx'G(T-t,x-x')0,Vx'), 

where 

is the advanced diffusion kernel which determines the propagator in the tree- 
graph solution here. 



4.3 Model D 



Model D has the same static action (23) as model C, but the kinetic coefficient 
is slightly different: 

( -AV^ " ^ 

a — 

-AiV 



V 



The matrix element an of this kinetic coefficient is the same as in model B. 



Therefore, the analysis of the previous section 4.2 is vahd here in full extent 
with the only replacement of the diffusion propagator G for the (f) field from 

(35) by the "hyperdiffusion" propagator G from relation (22) in the first equa- 
tion of set (34) and subsequently in the tree-graph expansion for (p^. Therefore, 
the constant is equal to Og. 



44 Model F 



In model F there are two complex conjugate fields: (pi = i/j, (p2 = 4'* ^"^^ ^ 
real field </73 = m. The static action is of the form: 

sf = I vv'r + fi^r + ^ + vMi^? ■ (36) 

D 2 
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The kinetic coefficient {a + (3) is determined by the matrices 



\^ 











-A^V^ 



(37) 



P 



ivs iv4ip 
—ivs —iv4ilj* 
—iv^ip* iv^ip 



(38) 



The static instanton is determined by the set of equations (25), (26) and (27) 
of the model C, where the filed is treated as a two-component one with the 
components ip*/^ 'ip/V^ and 0^ = ip*ip/2. Then the constant a = ap 
in relation (20) for this model was determined in section 4.2. The fluctuation 
determinant is also the same as in the model C and coincides with that of the 
massless static 0^ model. 

Equation (15) for the dynamic instanton is here 



'dt 

dm 
'dt 



(A^ + ivs) (^-VV + llV'PV' + V2mipj - iv^ip (m + V2\ip\^) , 
(A^ + ivs) (^-V^V'* + llV'TV'* + V2mip*^ + iv4ip* (m + V2\ip\'^) , 
-A^V^ (m + V2\'iIj\^) - iviV{i}*Vi} - i^Vi/j*) (39) 



with coefficients g = —1, V2 = — ix, V4 — * —iv^/u. The iterative solution can 
be written here in analogy with the model C using the advanced propagator 
G from (35). 



4.5 Model E 



Model E is F model with V2 — V3 — 0. Thus, it has been analyzed in the 

previous section. The constant in relation (20) is equal to ao due to V2 = 0. 
The static instanton for the field m is rrist = 0, but a nontrivial dynamic 
instanton for this field may be obtained from equations (39) due to mode 
coupling. 
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4.6 Model G 



In model G there are two real vector fields: let = 0a and </?3+a = rna^ where 
a — 1,2,3. The static action is of the form: 



and the kinetic coefficient is determined by the relations 

/ 



(40) 



a 





-A™V2 



(41) 



/3a6 = , (3aZ+b = V2eabc(f>c, Ps+aS+b = V2eabcmc , (42) 

where eahc is the completely antisymmetric tensor and all the indices a, b, c 
assume values 1,2,3. 

The static instanton sohition for this model is (j) = (pst, Tn = 0. Therefore 
g = —1 and the constant ac in relation (20) coincides with the static one for 
the massless 0*^ model: gq = oo- The fluctuation determinant is also the same 
as in the 0^ model. 

The dynamic-instanton equation (15) for this model can be solved as described 
in the preceding sections. 



4.7 Model H 



The static action of H model is of the form: 

Sh = I + f + ^ + ^ ^2m0^ + , (43) 

with the scalar fields (p and m and the transversal vector field v. The static 
instanton solution for the scalar fields turns out to be the the same as in 
model C, whereas Vgt — 0. Therefore, the fluctuation determinant has been 
also discussed in section 4.2. But in fact the critical dynamics is described by 
the simplified Hq model without the field m, whose static action is 

= I (V0)^ + + ^v^ . (44) 

This means that the constant an here as well as the fluctuation determinant 
are the same as in the usual 0^ model: 0^^ = 00. 
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The kinetic coefficient is given by the matrices 

/ \ V72 n \ / 



a 



-A^V2 
-A^V 



/3 



(45) 



V2V( 

y u -A^vy y-v2V(t) 

where the transverse projection operator for the vector field v is implied. 
Equation (15) for the dynamic instanton in the Hq model assumes the form 



d(l) 

'm 

dv 

'dt 



-AcV^v + i— (/.V^ 
u 



The iterative solution may be constructed with the use of the advanced prop- 
agators (22) and (35) [the propagator (35) for the vector field v taken in a 
convolution with the transversal projection operator, however]. 



5 Conclusion 



In this paper we have shown that the instanton method is applicable to large- 
order asymptotic analysis for all the near-equilibrium standard models B - 
H of critical dynamics in a form which is a generalization of the recently 
proposed dynamic- instanton approach in model A [5] . The factorial growth of 
the large-order coefficients of the perturbation expansion 

p[N] ^ ^, Ca^N' (46) 

has been proved for any quantity F calculable in the form of a perturbative 
series. The constants in expression (46) has been determined for all models 
mentioned. The exponent b in relation (46), however, depends on the quantity 
F considered. For any particular F it may be readily calculated with the use 
of the results of the present analysis. 

For instance, it may readily be proved that the exponent b in the e-expansion 
contribution to the dynamic index z is determined by the large-order asymp- 
tote of the expansion of the fixed-point value g^* and the second term of the 
perturbation series of the two-point correlation function. Therefore, in the 
0{n) symmetric dynamic theory the relation 



follows for the dynamic index z. 
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The fluctuation contribution to the amphtude factor C in relation (46) has 
been analyzed for all models B - H as well. 
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Appendix 



Let us outline the calculation of the fluctuation determinant at the instanton 
solution along the lines proposed for the model A [5] . Since we have extracted 
the u dependence from the dynamic action as a prefactor in Eq. (19), due to 
the stationarity conditions (12) - (14) the contribution from the fluctuations of 
u is independent of the fleld fluctuations and factorized as well. The remaining 
nontrivial fluctuation integral over 5ipa and (5</7^ may be written as 



Iexp< 6ip'^aab5^'^ 



OacWL + yOiab + Pab) T J h 7 



fD 



[aab + /^ab) 



5ipb5ipc5ipd 



^ r^ ^l^ab ^ 5^Pab 5S 

5(pc S(pb5(pd S(pc5(pd Sifb 



S^d}, (47) 



where f3^ is the fleld-independent part of the kinetic coefficient /3 (we remind 
that the kinetic coefficient a is assumed to be fleld-independent). The change 
of variables 



«6 V dt 



+ (Jbc-, — ? — dipd + -r — — dipd 



'S(pcS(pd 



5(pd Scpc 



(48) 



together with the antisymmetry of the kinetic coefficient f3, the instanton 
equation (15) and the second equality in relation (16) then allows to express 
the quadratic form of the exponential of the numerator of fluctuation integral 
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(47) as 
1 



hv,v'}^ ^ S^'^aabS^h 



dSipa 



dt 



+ 



{aab - Pab) 



S^S S(3ab SS 



1 5^S 

- o ^VstaJ 7 

2 Oipadcpb 



1 X 

(XPstb + o O^Oa^ 7 

2 5ipa5ipb 



5^ob 



(49) 



where Sipsta = 5(pa\t=T and 5(^oa = 5^a\t=to- To calculate the fluctuation 
integral pcrturbatively, we write — using the change of variables (48) — the 
normalization factor in terms of the perturbation expansion implied in (49) : 



d 

-^ac-KT + («a6 " Pab) ^bc 



dt 



^ SifistaKabSiPstb + ^ S(poaKabS(pob (50) 



with the same free-field part for the perturbative calculation as in (49). This 
representation gives rise to the propagators A12 = Aj^^ and the correlation 
function An determined by the operator Mq [defined by equation (7)], its 
transpose Mj and the correlation function of the Langevin source (2) as 



(51) 
(52) 



(A12),, {t, x; t', x') = [M,')^^ {t, x; t', x') , 
(An)„, (t, x; t\ x') = 2 (AuaA,,)^, (t, x; t\ x') , 



where A12 is the retarded Green's function of the operator Mq (and, corre- 
spondingly, A21 the advanced Green's function). 

Inspection of the Fourier transform of the propagator (Ai2)„b {t, x; t', x') shows 
that that in the time-space representation it is rapidly decaying at t — t' 00, 
which, together with the initial condition ^ lim^ {Ai2)^ {t, x; f, x') = SabS{x — 

x'), is basically all what is needed for an analysis similar to that for model A 
[5]. 

Let us introduce for the numerator of relation (47) the following decomposition 
of the integration variable 



(53) 



where 5ip consists of the standard Gaussian integration variable constructed 
with the aid of the advanced Green's function (A2i)^j and the solution h of 
the homogenous equation {Mq) hb{t,yi) = such that 5(^„(T, x) = 0, i.e. 



5Va = {An)abVb + {A2l)^^v'b + ha: 



(54) 
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where rja and r]'^ are fields vanishing at the boundaries (i.e. at |x| — > oo, i = T 
and t — to). This choice is accompanied by the prescription 



(55) 



for the auxihary field. 

The second term 5ip^ in the decomposition (53) is chosen as the solution of 
the homogeneous equation 



—OacW: + [<^ab — Pah) "? r r 



5(^J(i,x)=0 



(56) 



with the final value 5(f'^{T,:x.) = Sipstai'^) to conform to the boundary con- 
dition. The solution of equation (56) can be constructed perturbatively using 
the tree-graph expansion oi (po- 

The integration space for the normalization factor in the denominator of re- 
lation (47) is constructed with the aid of a decomposition similar to (53) 
with the same Sip and Sip'^. In this case, however, the second term Scp-^ in the 
decomposition may be expressed explicitly as the spatial convolution of the 
advanced Green's function (A2i)^5 (t, x; T, x') (note the fixed time argument) 
and the function 6ipsta{'^) and thus is the solution of the free-field counterpart 
of equation (56): 



-^ac^ + (ttafc - Pab) ^bc 



dt 



S^l{t,^) = 0. 



With the use of the asymptotic properties of the integration variables con- 
structed in this manner it may be shown that in the limit T — > oo the final 
(static) fiuctuations factorize completely (as in the case of model A): 



= y VS(fisteX.p (^~SipstaKabS(fistb 



X 



j VSifstl exp f —Sipsta 



6ipaSipb 



SiPstb 



dt 



X / / V6fV6^' exp{ 6ij'^aabS^h - 



VSifiVSi/j' exp I Si/j'^aabSi/j',, 
1 

+ ^^VOaKabS(Pob 
+ («a6 - Pab) 



dt 



Scpc 5(pb 



SipbSific 

^fob } > (57) 
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where 5({)q^ = SipJt=to- On the right-hand-side of (57) the first two factors 
yield the static fluctuation determinant including the degeneracy-lifting unit 
decomposition. A non-trivial dependence on the initial time instant to remains 
in (57) through the non-zero initial value of the fluctuation field 6(pQ. When 
the translation invariance with respect to time is restored by passing to the 
limit to —oo, however, the initial value of the fluctuation field vanishes: 
6(Pq —>■ due to the structure of the integration field (53) and the attenuation 
property of the advanced Green function. In a manner similar to that used 
for model A [5] it may then be shown perturbatively that in this limit the 
dynamic contribution to the fluctuation integral (57) tends to unity and thus 
we arrive at the fluctuation integral 



lim (5S 



[vSipst I exp f 6(psta J f 

J y 2 dLPadLpb 



completely determined by the static theory. 



-2 ^VstaKabSiPstb 



(5J 
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